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ABSTRACT 


A laser technique has been used to simulate processed 
output from a proposed spaceborne radio-astronomy 
interferometer experiment comprising two satellites 
separated by a controlled distance. The experiment 
effectively transmits the Fourier transform of the radio stars 
over the celestial sphere. The Fourier transform is a function 
of the distance between the two satellites and of the 
direction of the line connecting them. The existence of upper 
and lower bounds on the distance between the satellites 
creates voids in the Fourier transform space, causing 
resolution problems in the reconstruction of the star pattern 
on the ground from the Fourier transform. Because coherent 
optical techniques can obtain the two-dimensional Fourier 
transform, they can be used to simulate the proposed 
experiment. Optical filters can be made which represent the 
voids in the Fourier transform of the star pattern. An inverse 
transform shows the effects of the filtering, simulating the 
effects of the voids. This paper discusses optical filtering in 
general and the particular optical system used. Results are 
shown for 1-MHz and 10-MHz RF sources, representing the 
worst and the best resolution expected. 

Two conclusions are drawn. First, optical techniques are 
excellent for problems of this type. Second, because of the 
poor resolution in the 1-MHz range, 3,000 meters (the 
proposed maximum distance between the two satellites) is 
inadequate. 
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OPTICAL SIMULATION OF 

A RADIO-ASTRONOMY INTERFEROMETER EXPERIMENT 


by 

James P. Strong III 
Goddard Space Flight Center 


INTRODUCTION 


Coherent optical processing is a very powerful tool for making “spatial” frequency analyses of two- 
dimensional data. Optical techniques can produce a two-dimensional Fourier transform of the data which 
is a “spatial” frequency analysis. One can subsequently produce an inverse Fourier transform to reproduce 
the original data. Filters can be introduced which attenuate certain “spatial” frequencies of the Fourier 
transform. Taking the inverse Fourier transform will then show the effects of the filtering on the original 
data. This property of optical processing has made it useful in simulating a proposed spaceborne radio- 
astronomy interferometer experiment. 

This proposed experiment will attempt to make a 1- to 10-MHz RF map of the celestial sphere. For 
this purpose, two satellites will be placed side by side in the same orbit. Each will receive signals from 
astronomical radio sources, along with background noise. The received signal from one satellite (the 
maneuvering satellite) will be transmitted to the other (the reference satellite). The two signals, bandpass- 
filtered at the RF frequency under consideration, will be fed to a cross-correlator. Dr. E.J. Lueke, in a 
private communication, has derived the output of the cross correlator, C (N), as: 


C( N ) = JJJd( P) e~> (P ' N) dP 


( 1 ) 


P is a three-dimensional vector whose direction is from the two satellites to the RF sources (stars), 
as shov/n in Figure 1 . The magnitude of P is between zero and 1 . D( P) is the brightness function of the 
RF sources over a unit celestial sphere. D(P) has values only at |P| = 1 . N is a vector whose direction 
coincides with a line connecting the two satellites. Its magnitude is proportional to the distance between 
the two satellites and to the RF source' frequency under consideration. If one considers Equation 1, it 
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will be noted that C( N) is the three-dimensional 
spatial Fourier transform of the brightness func- 
tion D( P). When C(N) is relayed baek to earth, 
the inverse Fourier transform ean be taken to 
obtain: 


D( P) = 


There are several difficulties in obtaining the 
inverse Fourier transform. The one that concerns 
us the most is that N is limited in magnitude, 
since the two satellites can be neither infinitely far 
apart nor on top of one another. The proposed size 
of the antennas on each satellite limits the minimum separation to about 300 meters. It is expected that 
limited transmitter power for communications between the satellites will hold their maximum separation 
to about 3,000 meters. This limit places a restriction on the bandwidth of N, thereby reducing the 
resolution of Z)(P). 

An optical correlator is an ideal tool to exhibit the results of the limitations on |N|. It provides the 
experimenter with a photograph of a simulated RF celestial map showing the limits of resolution. 

The body of this paper will discuss the optically generated two-dimensional Fourier transform and 
how the transform can be used to simulate the magnitude limitations on N. It will then show, on a 
two-dimensional equivalent of an RF source map, the results of these limitations. 


/// 


C(N) e' (P ’ N > dN 


( 2 ) 



Figure 1 — Unit radius celestial sphere of 
interferometer experiment. 


OPTICAL FILTERING 


The optical experiment was set up as shown in Figure 2. where /, and/ 2 are the focal lengths of 
lens I and lens 2. respectively. If a transparency (for example, a star pattern) is placed at A. its 



Figure 2 — Coherent optical arrangement for obtaining Fourier transform and inverse. 
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two-dimensional Fourier transform will appear on a screen placed at B. If the screen is 
removed, lens 2 causes the inverse transform to appear on a screen placed at C. Thus, the Fourier trans- 
form of the Fourier transform at B produces the original image again at C. With this arrangement, one 
may place in the Fourier transform plane at B a mask which would attenuate certain regions of the 
transform plane. The results of this masking on the original image may be observed at C. This procedure 
is analogous to passing a signal through a filter and observing the output waveform. 

Let us now discuss the Fourier transform plane and see the results of masking. Figure 3 shows a 
two-dimensional Fourier transform plane. 


As in the more familiar one-dimensional 
case, the center (0,0) corresponds to zero fre- 
quency, or dc. A blank or clear transparency 
would produce a point of light at (0,0) in the 
Fourier transform plane, since the fundamental 
spatial frequency of such a transparency is zero. 
The higher the frequency components a trans- 
parency has, the farther out from (0,0) these com- 
ponents will plot in the Fourier transform plane. 
If an opaque mask with a circular hole of radius 
R were placed in the U-V plane and centered 
at (0,0), the result would be a “spatial low-pass 
filter” which would e liminate all frequencies in 
the U-V plane where Ju 2 + v 2 > R. To illustrate 
the transform plane and the filtering, consider 
Figures 4 and 5. Figure 4 is effectively a rotated 
square wave. Figure 5 is the Fourier transform 
of Figure 4, and in it can be noted the rings 
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Figure 3 — Fourier transform plane. 


corresponding to the fundamental frequency, the third, fifth, and seventh harmonics. The spot in 
the center corresponds to the average brightness or dc term. The rotated square wave has a Fourier 
series given by the equation 


A 0) = K 2 1 — cos (K x r) + -ycos (3K l r) — -~-cos (5^/-) + 


where K 2 and K x are constants and r is the distance from the center of the picture. The intensity of the 
brightness in a photograph such as Figure 4 is proportional to 


MO -)] 2 . 

If there is introduced at B a mask that eliminates the rings higher than the first harmonic (a low- 
pass filter), the result is 

A O') = K 2 [ 1 — cos (AT , r) ] 

and the intensity is 

Hr) = MOO] 2 = K 2 2 [ 1 — cos (AT, r) ] 2 
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Figure 4 — Rotated square wave. Figure 5 — Fourier transform of 

rotated square wave. 


Figure 6b shows this function compared to the original in Figure 6a. Thus we see that a removal of 
higher harmonics tends to blur sharp edges. If there is introduced at B a mask that eliminates the dc 
term and all rings higher than the third harmonic (a bandpass filter), the amplitude 


A (r) = K 2 


cos {K l r ) + 4-cos (3K x r) 


A photograph will record: 


I(r) = 


[ A O')] 2 = K 2 


-cos (Kjf) + j cos (3 K x r) 
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This is shown in Figure 7b. Note that areas which were dark in 7a (the original) are bright in 7b. When 
the dc component is removed, areas darker than the average brightness become negative in amplitude. 
When squared in the photographic process, the negative peaks become bright areas just as well as 
positive peaks. This is analogous to full-wave rectification. Thus “ringing” will be more likely to show 
up in photographs of the output than would be the case if an amplitude-sensitive device were used to 
record the output. 


The above examples show how coherent optical techniques can be used to study the effects of 
“spatial frequency filtering.” 


APPLICATION TO INTERFEROMETER EXPERIMENT 


Since, in Equation 1, C(N) represents the Fourier transform of D (P), the vector N represents 
points in the Fourier transform space. Limitations on N will therefore represent a type of “spatial 
frequency filtering.” The magnitude of N has been derived by Dr. E.J. Lueke as 


|N| =- f ± d - 
3x 10 8 


(3) 
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Figure 6 — Comparison of Figure 7 — Comparison of 

original square wave and original square wave and 

output through low - pass output through bandpass 

filter. filter. 

from which it can be seen that |N| is proportional both to the frequency / of the RF source and to the 
distance d between the two satellites. Thus for every RF source frequency there are a maximum |N| and a 
minimum |N| corresponding, respectively, to the 3,000-meter maximum distance and the 300-meter 
minimum distance: 

2irJ' 0 X 3,000 
3 X 10 8 

2irf 0 X 300 
3 X 10 8 


N 


N min 


These limitations on |N| correspond to masking the function C(N), whose magnitudes of |N| are 
unlimited, with a bandpass filter with minimum frequency N mjn and maximum frequency N max . 

To simulate this bandpass filter with the optical setup of Figure 2, the following experiment was 
performed. A two-dimensional equivalent of a unit radius celestial sphere was made. This consisted of a 
star pattern drawn within a circular area on a transparency. The resolution of the star pattern, or the 
maximum density of points in the pattern, corresponded to 20 stars per radian in the celestial sphere. 
The transparency was placed at A. Two masks were made in the configuration shown in Figure 8, an 
opaque disc being placed in the center to eliminate all spatial frequencies less than N min . One mask 
was made to simulate N max and N mjn for 10 -MHz RF sources and the other for 1 -MHz RF sources. 
The results, recorded on film at C (Figure 2), represent the best and worst resolutions obtainable from 
the proposed interferometer experiment. 
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Figure 8 — Mask for simulating two- 
dimensional bandpass filter. 



Figure 9 — Two-dimensional 
equivalent of unit radius 
celestial sphere. 


Figure 9 shows the original star pattern. Some 
short lines were placed into the picture to simulate 
sources such as gaseous clouds. Figure 10 shows 
the Fourier transform of the star pattern. The 
outer circle, I, and the middle circle, II, in Fig- 
ure 10 are N max and N min when the star pattern 
represents 10-MHz sources. The middle circle, II, 
and center circle, III, are N max and N min when the 
star pattern represents 1-MHz sources. Figures 1 1 
and 1 2 show the results of the filtering. Figure 1 1 
shows that bandpass limitations on 10-MHz RF 
sources have little effect on the resolution except 
for ringing around the line sources. This ringing 
may obliterate low-intensity point sources. Fig- 
ure 12 shows that the bandpass filter had severe 
effects on resolution of 1-MHz RF sources. Only 
the most intense point sources are discernible. The 
line sources are smeared in most areas of the 
picture. One might deduce the loss of resolution 
for 1-MHz sources by noting the small percent of 
the total Fourier transform information allowed 
to pass through between circles II and III in Fig- 
ure 10. 

To determine whether the resolution might 
be improved if frequencies less than N mjn were al- 
lowed, low-pass filter masks with N min = 0 were 
made (r, = 0 in Figure 8). Figures 13 and 14 
show the results. Figure 13 shows that, for 
10-MHz sources, frequencies less than N mjn have 
eliminated the ringing around the line sources. 
Figure 14 shows that, for 1-MHz sources, fre- 
quencies less than N mjn have little if any effect 
on the resolution. 


CONCLUSIONS 

The main conclusions concerning the 
proposed radio-astronomy interferometer 
experiment are the following: (1) Resolution of 
the reconstructed star pattern is degraded by 
having the lower limit on the satellite separation 
fixed at so high a value in the plotting of RF 
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sources frequencies near 10 MHz; (2) for resolution 
improvement in the 1-MHz range, the distance 
between the satellites must be much greater than 
3,000 meters. 

The simulation of this interferometer 
experiment is an excellent example of the ease 
with which coherent optics can be applied to a 
problem involving two-dimensional data. To do 
this problem on a digital computer would require 
long programs, large buffer memories, and input 
scanning devices. Coherent optics can do the job 
effectively with no program and no big buffers. 


Figure 10 — Fourier transform of star pattern 
of Figure 9, showing relative sizes of 
bandpass filters. 



Figure 1 1 — Reconstructed pattern 
of 10 -MHz RF sources through 
bandpass filter. 



Figure 12 — Reconstructed pattern 
of 1-MHz RF sources through 
bandpass filter. 
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Figure 13 — Reconstructed pattern 
of 10- MHz RF sources when 


Nmin 0. 


Figure 14 — Reconstructed pattern 
of 1 - MHz RF sources when 
N mjn =0. 




Appendix A 


Mathematical Treatment of Bandpass Restrictions 
on the Resolution of Point Sources 

The bandwidth problem mentioned in the introduction comes about from restrictions on the 
distance between the two satellites. The size of the antenna on the satellites limits the minimum distance 
to about 300 meters. Transmitter power limits the maximum distance to about 3,000 meters. 

The vector N directed parallel to the axis connecting the two satellites has a magnitude of 


|N| 


_ 27 r/ 0 d 
c 


(Al) 


which was derived by Dr. E.J. Lueke in a private communication. In this equation, f Q is the 
frequency of the RF source, c is the speed of light, and d is the distance in meters between 
the two satellites. Thus, with the proposed restrictions on satellite separation, |N| is restricted 
between the limits corresponding to 300 meters and 3,000 meters. 


In order to see mathematically the effects of these bandwidth limitations on the reconstruction of 
a point source, let us consider a one-dimensional case where the Fourier transform C(co) is given by: 


C(co) = 



<? ~i° JX 


dx 


(A2) 


In this equation, x is any variable, and to is the Fourier transform variable. Let the bandwidth limita- 
tion have the form shown in Figure A.1 , i.e., one having a ratio of minimum to maximum frequencies 
of 1/10, corresponding to the 300-to-3,000 meters distance limitation. 



Figure Al — Bandpass filter corresponding to N max and N m j n 
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Let / (x) be a Dirac delta function 5 (x — a) corresponding to a point source at x = a. Thus the Fourier 
transform of the point source is (from Equation A2) 


C(to) 


r + oo 

= / h(x — a 

J — OO 


) e /WJC dx = e /t 


The inverse Fourier transform is defined by 


f(*) 


-L 


4- oo 

C(co) e ioJX doj 


(A3) 


Taking into account the bandwidth restriction shown in Figure Al, the inverse transform for f(x) 
becomes 


f(x) 


-L 


+ b 

e 

b 


e iuJX doj 


/ 


+ b / 10 

e ~i wa 
- 6/10 


d cu 


which when solved becomes 


f(x) =2 b 


sin b (x — a) 

2b 

b (x — a) 
sm 10 

b{x - a) 

10 

(x - a) 


f (x) is plotted in Figure A2. From this figure it may be seen that the main brightness lobe has a width 

of about — . If another point source were at a distance less than -^ 7 - from a, the main lobes in 
b b 



Figure A2 — Reconstruction of f(x) = 6 (x-a) , showing 


results of bandpass filtering. 
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the reconstruction of f(x ) would overlap or blur together. This main-lobe concept can be expanded to 
the reconstructed star pattern over the celestial sphere. Consider again Figure A1 and let b = N max . 
Remembering from Equation A1 that 


the main-lobe width will be: 


N 


“T/o^max 

C 


27 T _ 277 - c = 1 0 5 

~T 2irf 0 d max 3,000 / 0 / 0 

c 

Thus the maximum resolution for a point source of 10 MHz will be 1/100 unit radius. Note that all 
dimensions are with respect to the radius of the unit celestial sphere. The 1/100 unit radius corresponds 
to a maximum of 100 stars per radian on the celestial sphere. For / 0 = 1 MHz, the maximum 
resolution will be 10 stars per radian. 


NASA-Langley, 1969 30 
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